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This presentation is a Moriond version of our recent paper (Kravtsov, Vikhlinin & Nagai'-') 
where we discussed X-ray proxies for the total cluster mass, including the spectral temperature 
(Tx), gas mass measured within rsoo (Mg), and the new proxy, Y^, which is a simple product 
of Tx and Mg. We use mock Chandra images constructed for a sample of clusters simulated 
with high resolution in the concordance ACDM cosmology. The simulated clusters exhibit 
tight correlations between the considered observables and total mass. The normalizations of 
the M500 — Tx, Mg — Tx, and M500 — relations agree to better than « 10 — 15% with the 
current observational measurements of these relations. Our results show that Y^ is the best 
mass proxy with a remarkably low scatter of only « 5 — 7% in M500 for a fixed Y,^, at both low 
and high redshifts and regardless of whether clusters are relaxed or not. In addition, we show 
that redshift evolution of the Y^ — M500 relation is close to the self-similar prediction, which 
makes Y^ a very attractive mass indicator for measurements of the cluster mass function from 
X-ray selected samples. 



1 Introduction 

The evolution of the cluster abundance is one of the most sensitive probes of cosmology. The 
potential and importance of this method have motivated efforts to construct several large surveys 
of high-redshift clusters during the next several years. However, in order to realize the full 
statistical power of the upcoming cluster surveys, it is paramount that the relation between 
cluster mass and observables and any potential biases are well known. 

Several cluster observables based on the galaxy velocities, optical light, X-ray observables 
such as luminosity, temperature, mass of the intracluster medium (ICM), and Sunyaev-Zel'dovich 
flux have been proposed as proxies of the total cluster mass (see a recent review by Voitl21). In 
this study we focus on the mass indicators derived from cluster X-ray observables. 

X-ray luminosity is the most straightforward mass indicator to measure observationally. 
However, is also the least accurate (internally) of all proposed X-ray proxies for Mtot with 
a large scatter EE and deviations of the slope of the — M relation from the self-similar 
prediction!^. The most common choice of mass proxy in the cluster cosmological studies has 
been the X-ray temperature of the intracluster gas ' . The scatter in the M — T^ relation is 
smaller compared to that in the L ^ — M relation (the upper limit from observations is ~ 15% in 
M for fixed T for relaxed clusters^!). In general, existence of a tight relation such as M — Tx 
indicates that clusters are regular population of objects with their global properties tightly 
related to total mass, and scatter caused bv secondary effects such as substructure in the ICM, 
non-gravitational processes, and mergers'^. More recently, gas mass was used as a proxy for 



Mtot ' '■ The practical advantage of Mg over Tx is that it can be measured from the X-ray 
imaging alone. Also, Mg can be expected to be less sensitive to mergers which should translate 
into smaller scatter in the Mg — M relation. The caveat is that trend of gas mass with cluster 
mass and evolution with redshift are not yet fully understood. 

The use of clusters as efficient probes for precision cosmology puts stringent requirements on 
observable cluster mass proxies: 1) tight, low-scatter correlation between the proxy and mass, 
with the scatter insensitive to mergers etc., and 2) simple power-law relation and evolution which 
can be described by a small number of parameters and be as close as possible to the prediction 
of the self-similar model. The last point is cmcia l to ensure that the self-calibration strategies 
for analyses of large cluster surveys ^" 4 | 15 | 16 | 17 | 18 | successful. This is because self-calibration 
is powerful when cluster scaling relations and their evolution have a simple form which can be 
parameterized with a small number of parameters. 

In general, a mass proxy does not have to be a single cluster property, such as L^, or Mg. 
Any physically-motivated combination of these variables that is expected to be tightly related 
to cluster mass can be used to construct a valid mass indicator. A hint for a better X-ray mass 
proxy is provided by recent studies based on cosmological simulations of cluster formation 
which show that integrated SZ flux, Y^z is a good, robust mass indicator with low scatter in 
the Ysz — M relation, regardless of the dynamical state of the cluste r. In a ddition, the Isz — M 
relation exhibits a simple, nearly self-similar evolution with redshift^^^^. The physical reason 
for the robustness of the SZ flux is straightforward: Ysz is directly related to the total thermal 
energy of the ICM and thus to the depth of the cluster potential well. 

Here we show that a similar robust, low-scatter mass indicator can be constructed using X- 
ray observables. The indicator, which is simply the product of the X-ray derived gas mass and 
average temperature, 1^ = Mg Tx, correlates strongly with cluster mass with only ~ 5 — 8% 
intrinsic scatter. The 1^ — Mtot relation is robust to mergers, in the sense that even for disturbed 
unrelaxed systems it gives unbiased estimates of mass with the statistical uncertainty similar to 
that for relaxed systems. In addition, we show that evolution of the slope and normalization 
of the Yx — M relation is nearly self-similar. These properties make Yx particularly useful for 
measurements of cluster mass function using X-ray surveys. 



2 Mass Proxies 



Physical properties of virialized systems, such as clusters, are expected to correlate with their 
total mass. For example, in the self-similar model^^ the cluster gas mass is expected to be 
simply proportional to the total mass, Ma^ = CugMg^/^^, where masses are determined within 
a radius enclosing a certain overdensity Ac with respect to the critical density of the universe 
at the epoch of observation, pcrit{z), and CMg is a constant independent of cluster mass and 
redshift. The self-similar relation between cluster mass and temperature is E(z) Ma^ = CtT^^'^- 
Here the function E{z) = H{z)/Hq for a flat cosmology with the cosmological constant assumed 
throughout this study is given by E{z) = ($7a/(1 + z)^ + ^^a)^''^- 

The SZ flux integrated within a certain radius, Y^z, is proportional to the total thermal 
energy of the ICM gas and thus to the overall cluster potential, which makes it relatively 
insensitive to the details of the ICM physics and merging. Isz is proportional to the ICM mass 
and gas mass-weighted mean temperature, Y^z oc Mg^Ac^m- The self-similar prediction for the 
Ysz — M relation is 

E{zf/'MA^ = Cy,,Yli\ (1) 

Cosmological simulations show that Ysz is a good, low-scatter cluster mass proxv and that 
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Ysz — M relation form and evolution are close to the self-similar prediction ^ozu_ Qiven 
the good qualities of Ysz as a mass proxy, it is interesting whether a similar indicator can be 




constructed from the X-ray observables, which could be used in studies of the X-ray cluster 
abundances. The simplest X-ray analog of Ysz is 

i; = Mg,A.T., (2) 

where Mg^Ac is the gas mass derived from the X-ray imaging data (it is measured within a 
radius enclosing overdensity Ac), and Tx is the mean X-ray spectral temperature. is measured 
excluding the central cluster region, which can be achieved with moderate angular resolution 
X-ray telescopes (< 15" FWHM). To excise the central regions is desirable because the observed 
cluster temperature profiles show a greater degree of similarity outside the core^^, and also 
because this makes the spectral temperature closer to the gas mass averaged T^. 



3 Mock Chandra Images and Analyses of Simulated Clusters 

To test the quality of Tx, Mg, and Yx as the total mass proxies, we take the following approach. 
We use the output of the high-resolution cosmological simulations of clusters in a wide mass 
range to predict their X-ray emission maps. These maps are convolved with the response of the 
Chandra telescope to create realistic mock "observations" of these clusters. The mock data are 
used as an input to the actual X-ray data analysis pipeline and measure the proxies (Tx, Mg, 
and Yx) as they would be derived by observers. 

The simulations and analysis procedure is fully described in Kravtsov et al. The cosmo- 
logical simulations accurately follow the cluster gro wth from inital conditions using Adaptive 
Refinement Tree (ART) A^-body-|-gasdynamics code^^, a Eulerian code that reaches the high 
dynamic range required to resolve cores of mass halos. The peak formal resolution achieved 
in these simulations, ~ 3.66 kpc, is sufficient to resolve halos of individual galaxies and to 
follow not only dissipationless dynamics of dark matter and gasdynamics of the ICM, but also 
star formation, metal enrichment and thermal feedback due to the supernovae type II and type 
la, self-consistent advection of metals, metallicity dependent radiative cooling and UV heating 
due to cosmological ionizing background. 

Generation and analysis of the mock X-ray data reproduces all the effects associated with 
the Chandra response and reconstruction of the 3D ICM properties from the observed projected 
X-ray spectra. The only exception is that we ignore complications present in reduction of the 
real Chandra data related to background subtraction and spatial variations of the effective area 
(i.e., we assume that accurate corrections for these effects can be applied to the real data and 
any associated uncertainties are included in the reported measurement errors). 
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Figure 2: Relation between the X-ray spectral tem- 
perature, Tx, and total mass, Msoo- Separate symbols 
indicate relaxed and unrelaxed clusters, and also z = Q 
and 2 = 0.6 samples. The dashed line shows the power 
law relation with the self-similar slope fit to the entire 
sample, and the dotted lines indicate 20% scatter. 
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Figure 3: Correlation between gas mass and total mass 
of the clusters. Both masses are measured within rsoo. 
The meaning of the symbols and lines is the same as 
in Fig. 121 The dotted lines indicate 15% scatter. 



Our simulations include clusters in a wide range of mass and dynamical state. We use 
simulation outputs at redshifts z = and z = 0.6 to test the evolution in the mass vs. proxy 
relations. Examples of the mock Chandra images are shown in Fig^ 

4 Comparison of Mass Indicators 

Figure HI shows that the slope and evolution of the M500 — relatiorf are quite close to the 
self-similar model. There is a ~ 20% scatter in M500 around the mean relation and much of 
the scatter is due to unrelaxed clusters. Note also that the normalizations of the M500 — 
relation for relaxed and unrelaxed systems are somewhat different: unrelaxed clusters have 
lower temperatures for a given mass. This may seem counter- intuitive at first, given that one 
can expect that shocks can boost the ICM temperature during mergers. However, in practice 
the effect of shocks is relatively small^^. The main source of the bias is that during advanced 
mergers the mass of the system already increased but only a fr actio n of the kinetic energy of 
merging systems is converted into the thermal energy of the ICmI^^. 

The M500 — Mg relation (Fig. |3J has a somewhat smaller scatter (~ 10 — 12%) around the 
best fit power law relation than the M500 — Ty^, but its slope is significantly different from the 
self-similar prediction — we find M500 oc Mg-^^~^-^'^ compared to the expected M500 oc Mg. This 
is due to the trend of gas fraction with cluster mass, /gas = Afg/Msoo oc M^qq""'^ present for 
both the simulated clusters in our sample!^ and for the observed cluster^^. The normalization 
of the Msoo — Mg relation evolves only weakly between z = 0.6 and z = (yet, the evolution is 
statistically significant and it reflects slow evolution of the gas fraction with time^^^ . 

The M500 — Yx relation (Fig. 2)) has the smallest scatter of only ~ 5 — 7%. Note that this 
value of scatter includes clusters at both low and high-redshifts and both relaxed and unrelaxed 
systems. In fact, the scatter in M500 — Yx for relaxed and unrelaxed systems is indistinguishable 
within the errors. Note also that the figures include points corresponding to the three projections 



"Hereafter, M500 is the total mass within the sphere corresponding to the mean overdensity of 500 relative to 
the critical density at the cluster redshift. 
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Figure 4: 1^ — M500 correlation. The meaning of the 
symbols and lines is the same as in Fig.H The dotted 
lines indicate 8% scatter. 
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Figure 5: Fractional deviations in Tx and Mg for fixed 
Mtot from their self-similar relations. Solid and open 
circles show clusters at 2: = and z = 0.6, respec- 
tively. The deviations for Mg and Tx are generally 
anti-correlated. A similar anti-correlation exists in the 
trend with redshift. 



of each cluster. Figure 0] shows that the dispersion in the projected values of for each given 
cluster is very small, which means that 1^ is not very sensitive to the asphericity of clusters. 
Remarkably, the scatter of the M500 — relation, which involves direct X-ray observables, is as 
small as that in the M500 — i^sz relation 7% for our sample). 

The comparison of the mass proxies, clearly shows that 1^, the product of gas mass and X-ray 
spectral temperature, is more robust and self-similar mass indicator than either of these X-ray 
observables. Why is the product better than its parts? The answer is obvious from Figure [5] 
where we show the residuals of temperature and gas mass from their respective relations with 
total mass. Figure^lshows that the clusters with temperatures lower than the mean temperature 
for a given total mass tend to have gas mass higher than the mean, and vice versa. Note also 
that there is some redshift evolution between z = and z = 0.6 — more clusters have negative 
deviations of temperature and positive deviations of measured gas mass at z = 0.6 compared to 
z = 0. This redshift evolution is thus in the opposite direction for the gas mass and temperature 
deviation. The measured Mg systematically increases at higher z for a fixed total mass because 
high-z clusters are less relaxed on average. For unrelaxed clusters, the ICM d ensi tv distribution 
is non-uniform which results in overestimation of Mg from the X-ray data^^. Some of the 
decrease of Mg at lower z may be due to continuing cooling of the ICM which decreases the 
mass of hot. X-ray emitting gas. The anti-correlation of residuals and opposite evolution with 
redshift for gas mass and temperature is the reason why the behavior of their product, on 
average, has smaller scatter and is closer to the self-similar expectation in both the slope and 
evolution. 



5 Discussion and Conclusions 

We presented comparison of several X-ray proxies for the cluster mass — the spectral tempera- 
ture Tx and gas mass Mg derived from the X-ray data within rsoo, and the new proxy, Y^, defined 
as a simple product of and Mg. Analogously to the integrated Sunyaev-Zel'dovich flux, 1^ is 
related to the total thermal energy of the ICM. To test these mass proxies, we use mock Chandra 
"observations" of a sample of clusters simulated in the concordance ACDM cosmology. 



The main result of this study is that 1^ is a robust mass indicator with remarkably low 
scatter of only ~ 5 — 7% in M^qq for fixed 1^, regardless of whether the clusters are relaxed 
or not. In addition, the redshift evolution of the 1^ — M500 relation is close to the self-similar 
prediction given by equation which makes this indicator a very attractive observable for 
studies of cluster mass function with the X-ray selected samples. 

The Tx — M500 relation has the largest scatter (~ 20%), most of which is due to unrelaxed 
clusters. The unrelaxed clusters have temperatures biased low for a given mass because a certain 
fraction of the kinetic energy of merging systems is still in the form of bulk motions of the ICM. 
The Mg — M500 relation shows an intermediate level of scatter, ~ 10 — 12%. This relation does 
not appear to be sensitive to mergers. It does, however, exhibit significant deviations from self- 
similarity in its slope, which is due to the dependence of gas frac tion within rsoo on the cluster 
mass^^ (a similar dependence exists for the observed clusters^^). 

Generally, all the observable-mass relations we tested demonstrate a remarkable degree of 
regularity of galaxy clusters as a population. T^, Mg, and all exhibit correlations with M500 
which are close to the expectation of the self-similar model, both in their slope and evolution 
with time, within the uncertainties provided by our sample. The only exception is the slope of 
the Mg — M500 relation. 

Given that our analysis relies on cosmological simulations, it is reasonable to ask whether 
the simulated clusters are realistic. Although simulations certainly do not reproduce all of the 
observed properties of clusters, especially in their core regions, the ICM properties outside the 
core in simulations and observations agree quite well. We illustra te t his in Fig. 1^1 which shows 
that the Mg — Tx relations for simulated and observed clusters Clearly, both simulated 
and observed clusters exhibit tight correlations between Mg and which agree remarkably in 
their slope {Mg oc T^' ''^) and normalization. The normalizations derived from simulated and real 
clusters agree to ~ 10%, while slopes are indistiguishable and both deviate significantly from the 
expected self-similar value of 1.5. This is a consequence of si gnif icant trends in the ga s fra ction 
with cluster mass, Mg/Msoo oc M^oo^°-2^ for both simulated!^ and observed clusters'!^. The 
deviations from the self-similar model also manifest themselves in the absence of any noticeable 
evolution with redshift?. Interestingly, the real clusters show a similarly weak evolution in the 
Mg - Tx relation^. Fi gure El shows that the likely explanation is that the clusters at z = 0.6 
tend to be colder for the fixed Mtot but have higher estimated Mg then their counterparts at 
2 = because they are less relaxed. 

A similar level of agreement between the simulations and latest Chandra measurements ex- 
ists also for the total mass vs. temperature relation, Msnn — T^. In fact, the normalization for 
our simulated sample agrees with the observational results'^ to ~ 10%. This is a considerable 
improvement over the situation of just several years ago when there was ~ 30 — 50% discrepancy 
between observational measurements and cosmological simulations The M — normal- 

ization was revised both in simulations and observations due to (1) inclusion of more realistic 
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physics in cosmological simulations (e.g., radiative cooling and star formatio n^ I 'J^ i Q '^ i (2) im- 
proved analyses of observed clusters using more realistic gas density profiles^ ^ ^ (3) more 
reliable measurements of the cluster temperature profiles and (4) the use of uniform 
definition of Tx in observations and in simulations analyses^^^^^^^^. The remaining systematic 
10% differen ce observ ed at present is likely caused by non-thermal pressure support from bulk 
gas motions I 39 | 40 | 4l ] ig unaccounted for by the X-ray hydrostatic mass estimates. In 

Figure [3 we compare the — M500 relation for the simulated clusters and for the Chandra 
sample^. The observed clusters show a tight correlation with the slope close to the self-similar 
value. There is ~ 15% difference in normalization, likely explained also by neglecting the tur- 
bulent pressure support in the Chandra hydrostatic mass estimates. The excellent agreement of 



''Note that Mg in Fig|^is not multiplied by the E{z) factor unlike the total mass in Figl^andlH 
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Figure 6: Relation between X-ray spectral temperature 
and gas mass for the relaxed subsample of simulated 
clusters (circles) and for a sample of relaxed Chandra 
clusters (stars). Both gas mass and temperature are 
the quanities derived from analysis of real and mock X- 
ray data. The error bars in the Chandra measurements 
are comparable to the symbol size and are not shown 
for clarity. The dashed line shows the best fit power 
law relation with the slope 1.75. 



Figure 7: Kc — M500 relation for the 2 = sample 
of the simulated clusters (circles) and for a sample of 
relaxed Chandra clusters. The gas masses for the sim- 
ulated clusters are appropriately rescaled (see caption 
to FiglSJ. The dot-dashed line shows the best fit power 
law relation for the simulated clusters with the slope 
fixed to the self-similar value of 3/5. The dashed line 
shows the same best fit power law, but with the nor- 
malization scaled down by 15%. 



simulations and observations in terms of the relation between the two X-ray observables used to 
compute (Mg — T^) and a relatively good agreement in the — M500 and 1^ — M500 relations, 
gives us confidence that the results presented here are sufficiently realistic. 

Our results show that is clearly most robust and most self-similar X-ray cluster mass 
indicator. The biases existing in mass estimates based on Mg and anti-correlate both for 
a given redshift and in terms of evolutionary trends (see Figure |SJ . This explains why their 
product, Yx, is a better mass indicator than Tx and Mg individually. The quality of Y^ compares 
well to that for the actual three-dimension integral of the ICM thermal energy (proportional to 
Isz) in terms of its low scatter and self-similarity. Y^ may prove to be an even better mass proxy 
than igz, given that we use ideal 3D measurement of the latter while reproducing the actual 
data analysis for the former. Note also that Y^z is more sensitive to the outskirts of clusters, 
because it involves gas mass- weighted temperature (as opposed to the spectral temperature more 
sensitive to the inner regions), and thus should be more prone to projection effects. 

Note that Y^ is also an attractive mass proxy from the data analysis point of view. First, it 
reduces observational statistical noise by combining the two independently measured quantities, 
Mg and Tx, into a single quantity. For example, a 10% measurement uncertainty in Tx translates 
into a ~ 15% mass uncertainty through the M — Tx relation and only 6% uncertainty through 
the Y^ — M relation. 1^ is also less sensitive to any errors in the absolute calibration of the 
X-ray telescope because the biases in the derived Tx and Mg tend to anticorrelate. 

The robustness and low scatter make Y^ an excellent mass indicator for observational mea- 
surements of cluster mass function at both z = and higher redshifts. The necessary data — 
an X-ray brightness profile and a wide-beam spectrum excluding the core — are easily obtained 
with sufficiently deep observations with Chandra, XMM-Newton, and Suzaku (for low-redshift 
clusters). The small scatter and simple, nearly self-similar evolution of the 1^ — M relation hold 
promise for the self-calibration strategies for future large X-ray cluster surveys. 
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